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Abstract 

Qh' In conventional general relativity, the expansion rate H of a Robertson- Walker 

Q ■ universe is related to the energy density by the Priedmann equation. Aside from the 

present day, the only epoch at which we can constrain the expansion history in a 
model-independent way is during Big-Bang Nucleosynthesis (BBN). We consider a 
simple two-parameter characterization of the behavior of H during BBN and derive 
constraints on this parameter space, finding that the allowed region of parameter 
r> I space is essentially one-dimensional. We also study the effects of a large neutrino 

' asymmetry within this framework. Our results provide a simple way to compare an 

alternative cosmology to the observational requirement of matching the primordial 
abundances of the light elements. 



1 Introduction 

Modern cosmology boasts a best-fit model which is in good agreement with a variety of 
data. This model features a homogeneous and isotropic spatially-flat universe comprised 
of approximately 5% baryons, 25% cold dark matter, and 70% vacuum energy, with a 
radiation energy density of about 10~^ at a temperature 2.7° K. Along with a nearly scale- 
free spectrum of adiabatic density perturbations, these ingredients provide an accurate 
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match to observations of the recent expansion history large-scale structure P], the 
cosmic microwave background (CMB) 0, and the abundances of light elements as predicted 
by Big-Bang Nucleosynthesis (BBN) [||. 

As successful as the best-fit model has been, it falls short of achieving an aura of 
inevitability due to certain problems of naturalness. Foremost among these are the two 
cosmological constant problems: why the vacuum energy is so small in comparison to its 
expected value, and why its precise magnitude is so close to that of the matter energy density 
today 0. In addition to these unsolved problems, there are problems for which a popular 
solution exists but is lacking a firm experimental footing: the flatness and horizon problems, 
which may be solved by inflation, and of course the nature of the dark-matter itself. Finally, 
amidst the large successes there are small failures of the best-fit model in fitting the data, 
especially in the detailed matching of structure on small scales to the predictions of cold 
dark matter [Q. Meanwhile, numerous alternatives to conventional cosmology have been 
studied. To take one recent example, brane-world models with large extra dimensions can 
give rise to a variety of departures from four-dimensional general relativity, with important 
consequences for the early universe 0, ^ . 

Given this situation, it is worth examining as carefully as possible the theoretical as- 
sumptions of the standard cosmological model, if for no other reason than to reassure 
ourselves that its apparent fine-tunings are not pointing toward a radically different under- 
lying framework. In this paper we consider the empirical evidence relating to the Fried- 
mann equation, the dynamical relation in general relativity between the expansion rate of 
the Universe and the energy density in it. 

The fundamental quantity in general relativity is the space-time metric. We have some 
empirical evidence for the form of the metric from astronomical observations. We know 
that the Universe is expanding 0, i.e., the physical volume of the Universe is increasing 
with time. We also know that the cosmic microwave background is isotropic to a part in 
10^ [|TU|. Therefore, it seems reasonable to assume that the spatial metric is isotropic. If we 



further assume that all observers (matter) in this expanding isotropic Universe see the same 
CMB isotropy, then it follows |TT| that the Universe is homogeneous (on scales much larger 
than that of the largest collapsed structure). The metric is then uniquely determined. The 
form of this metric, called the Robertson- Walker (RW) metric, is 



ds^ = -df + a^(t) 



(1) 



- kr 

where a{t) is the scale factor and k G { — 1,0,+!} is the curvature parameter. The Fried- 
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mann equation (derived from general relativity for the RW metric) is 



,,2 8ttG k 



a J 



(2) 



H = a /a is the Hubble expansion rate and p the energy density of the universe. In fact 
d^) is only the 00 component of Einstein's equations applied to (0); we also have another 
independent gravitational equation, 
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as well as an equation of energy-momentum conservation, 



-3{p + p)- 



(3) 
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where p is the pressure. In seeking to test this framework, we might hope to seek consistency 
relations among these equations, which could be compared with data. Unfortunately, it 
is always possible to find an energy density p{t) and pressure p{t) which satisfy 
Specifically, any given expansion history a{t) and curvature k is compatible with these 
equations if we choose 
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The crucial point here is that a perfectly smooth component of energy and pressure cannot 
be detected in any way except for its influence on the expansion rate, so we have no 
independent constraint on such a source. While invoking a perfectly smooth component to 
fit an arbitrary behavior of the scale factor might seem like cheating, this is exactly what we 
must do to reconcile the strong evidence in favor of spatial flatness with the similarly strong 
evidence that the amount of clustered matter falls far short of the critical density in the 
current universe. In this sense, the set of equations (0-^) are, strictly speaking, untestable, 
without some prior expectation for the nature of p and p. 

Instead, we can characterize what observations can tell us about the behavior of the 
scale factor in a model-independent way, so that any specific alternative theory can be 
straightforwardly compared with the data. There are only two eras of the universe's history 
in which such a characterization is possible: the recent universe, and the BBN era. In the 
recent universe, information about the behavior of the scale factor can be obtained in a 
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variety of ways, most directly by comparing the distance of faraway objects to their redshifts, 
as is done in the supernova studies that first revealed the acceleration of the universe 
During BBN, the expansion rate directly affects the relative abundances of light elements 
produced (as reviewed below), so that the primordial abundances are a powerful constraint 
on alternative expansion histories. 

Other important observable phenomena which are affected by the expansion rate do not 
offer such a direct test, since they typically involve the local behavior of gravity (evolution 
of perturbations) in addition to its cosmological behavior (evolution of the scale factor). 
Examples include the growth of large-scale structure and the imprinting of temperature 
anisotropics on the CMB. A theory which predicts deviations from the Friedmann equation 
could generically predict deviations from the predictions of general relativity on other scales. 
Consequently, it is difficult to use the growth of structure and CMB anisotropies as model- 
independent tests of the Friedmann equation, although any fully-specified alternative theory 
might be tightly constrained by these phenomena. However, if one assumes that the local 
behavior of gravity as predicted by general relativity is correct, then structure formation 
and CMB anisotropies are a powerful probe of the expansion history (and hence dark energy 
evolution) after the last scattering epoch [|12 . 



In this paper we will study what kinds of expansion histories in the early universe are 
compatible with the BBN explanation of the light element abundances. Other works have 
constrained the energy density [|l^ or value of Newton's constant ||14| during BBN, assuming 
the Friedmann equation, or have derived BBN constraints on specific scalar-tensor theories 



15, 16 1 or have put limits on alternative cosmologies under the assumption that the universe 



has undergone a consistent power-law evolution from very early times fl^. An estimate of 
the constraints imposed by BBN (by considering the change in the neutron-proton freeze- 
out temperature) and structure formation on cosmologies where the energy density in the 
universe varies as some power of H has also been calculated |p. Our approach will be not to 
assume any specific behavior of the scale factor over long periods, but instead to introduce a 
two-parameter family of evolution histories which we take to be valid only in the vicinity of 
BBN. We find that a variety of alternative cosmologies can be consistent with observations, 
although they comprise essentially a one-dimensional region in our two-dimensional space 
of possibilities. 
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2 Parameterizing the scale factor during BBN 



We will consider theories in which the field equations for the metric may be different from 
those of general relativity, but we assume that test particles will still follow geodesies of this 
metric. From this assumption it follows that the energy of a relativistic particle redshifts 
as 1/a. The photon temperature will be proportional to 1/a in the absence of entropy 
creation; however, annihilation can act as a significant entropy source. We therefore use 
the neutrino temperature T as a measure of the scale factor, since T oc 1/a to excellent 
accuracy during the epoch under consideration. Of course these statements rely on our 
decision to only consider changes in the gravitational dynamics, not any particle-physics 
processes. 

The process of conventional BBN spans temperatures between the freeze-out of weak 
interactions, Tf = 1 MeV, and the synthesis of helium, Tbbn = 60 keV; this corresponds to 
a change in the scale factor by slightly more than one order of magnitude. It will therefore 
be reasonable to approximate the expansion rate during this interval, so long as it is not 
wildly oscillating or somehow finely tuned, as a simple power law: 

Any expansion history is parameterized by the coefficient Hi and the exponent a. (One 
can view this as a Taylor expansion in \og{H), to first order in log(T).) We do not extend 
this parameterization beyond the BBN epoch, nor are we proposing any models in which 
(0) is predicted (but note that all the models in |^, are well described by this expansion 
law); this is a phenomenological study of the expansion rate during BBN. It should also 
be noted that there might exist models which are not accurately described by power-law 



expansion during the BBN era in such cases it would be necessary to examine each 
model individually. 

In standard BBN weak interactions freeze out close to T/ = 1 MeV, at which point 
n/p{Tf) ^1/6. At temperatures lower than Tf, the only change in the neutron to proton 
ratio occurs due to free decay of neutrons with a lifetime of = 887 sec. At temperatures 
much lower than freeze-out, it becomes possible to synthesize helium in appreciable quan- 
tities. This happens at a temperature of Tbbn = 60 keV, and essentially all the remaining 
neutrons (one per seven protons) are bound into helium nuclei. It should be noted that 
this temperature (Tbbn) is much lower than the binding energies of all the light elements. 
Clearly, the large photon to baryon ratio (which implies a large photo-dissociating back- 
ground) is partly responsible for this |jl8[. The rate at which n —>■ ^He conversion occurs 
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depends on the destruction rate of deuterium which in turn depends on the binding energy 
of deuterium, the formation rate of tritium and the baryon to photon ratio. In standard 
BBN, it turns out that the conditions for hehum synthesis are just right when the hehum 



formation rate is larger than 1/rn ||T9l. This picture breaks down in the unconventional 
expansion histories we consider. 

An understanding of how alternative models affect light-element abundances can be 
expressed in terms of the behavior of three characteristic quantities: the freeze-out temper- 
ature Tf (which sets the initial neutron/proton ratio), the time interval between freeze-out 
and nucleosynthesis, A^bbn = ^bbn — tf (during which neutrons decay), and the Hubble 
parameter at nucleosynthesis, -ffBBN (which affects the efficiency with which neutrons are 
converted to helium, and thus determines the deuterium and lithium abundances). Note 
that AtBBN ~ ^BBN, since tf « teBN; this will continue to be true in the non-standard 
models we study. Assuming only that H is monotonically decreasing even prior to the 
epoch when this formula is relevant, we will have 

1 



^BBN 



aH] 
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The predictions of Big-Bang nucleosynthesis are given as ratios of abundances of the 
light elements - neutrons (n), protons (p or H), deuterium (D), '^H, ^He, ^He and ^Li. One 
quotes the ratio of the number density of the light elements to the number of protons, 
except in the case of helium where the accepted practice is to quote its mass fraction (Ip). 
The absolute density of these light elements is set by the absolute density of baryons. The 
only input parameter in standard BBN is the baryon density quoted as its ratio [r]) to the 
number density of photons. It is useful to define t]io = lO^^t] since the number density 
of baryons (compared to that of photons) is very small. Except for a brief period when 
electrons and positrons annihilate, r] remains constant (provided no other entropy changing 
interactions come into equilibrium). 

We will adopt observationally allowed ranges of light elements ratios in keeping with 



the inferences (but slightly on the conservative side) of Olive et al. pO[. The ranges are 
0.228 <Yp < 0.248, 2 < 10^ x D/H < 5, and 1 < 10^° x ^Li/H < 3. The low deuterium 
value has been adopted since it seems to be favored by current data pT| . 

The baryon density is not as well constrained observationally as the ratios of light 
element abundances. One reason for this is that we can only inventory the baryons that 
we see and most of the baryons could be dark. Even putting a lower bound on the baryon 
density is a difficult task. Persic and Salucci P2| estimated the density of baryons in the 
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Universe (in terms of r/io) to be r/io > 0.3. In converting the baryon density to r] we have 
assumed that rj has remained constant since BBN (as it does in the standard model), and 
that -f^o > 60 km s^^ p3| where Hq is the present expansion rate of the Universe. Fukugita 
et al. also did a baryon inventory and their bound on the baryon density translates to 
^710 ^ 1-8) a much higher lower bound on the baryon density. We will use both these lower 
bounds as a way of showing how the allowed region changes with the bound on the baryon 
density. 

Without recourse to a specific theory of gravity it is not possible to put an upper bound 
on the baryon density. For our purposes, we will assume the following ranges for the 
baryon density - a tight bound corresponding to 1 < r/ig < 10 and a broader range with 
0.5 < ?7io < 50. We will also comment on what happens if the upper bound on the baryon 
density is increased. 



3 Helium constraints 

The requirement that about 24% helium by mass be synthesized imposes severe constraints 
on the expansion history for temperatures around Tbbn- To understand these constraints, 
it proves useful to look at the effects of changes in the parameters Hi and a of (0) in the 
vicinity of the standard model. We first concentrate on the effect of changing Hi at fixed 
a. As mentioned, the important considerations are the freeze-out temperature Tj, and 
the time interval between freeze-out and nucleosynthesis A^bbn ~ ^bbn- (In the vicinity 
of standard BBN, small changes in the expansion rate -ffBBN will not be important, since 
essentially all neutrons are converted to helium.) Increasing Hi increases the expansion 
rate at every temperature; thus, freeze-out will occur earlier (Ty will be higher) leading to 
a larger initial neutron/proton ratio. At the same time, the faster expansion rate leads to a 
decrease in ^bbn, leaving less time for neutrons to decay. Both these effects go in the same 
direction and hence the neutron/proton ratio at Tbbn increases with increasing Hi, leading 
to a higher helium abundance. 

One can analyze the effect of changing a at fixed Hi in a similar manner. Again, the 
relevant considerations in the vicinity of standard BBN are the behaviors of Tj and A^bbn- 
Since Hi corresponds to a temperature close to freeze-out and before helium synthesis, in- 
creasing a means that the expansion rate will be lower during nucleosynthesis. In particular 
this implies that -?/bbn will be lower (i.e., helium synthesis happens later), and hence that 
the time interval AtBBN oc I/T/bbn will be longer. Consequently, increasing a at fixed Hi 
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gives neutrons more time to decay between Tf and Tbbn, working to decrease the helium 
abundance for values close to the standard ones. 

We can therefore balance the effects of increasing Hi against those of increasing a, 
to obtain a constant helium abundance for small deviations from the standard picture. 
In Figure |I| we plot contours of constant ^He in the a-Hi plane. This figure reveals the 
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Figure 1: Contours representing constant helium in the a-Hi plane, for Yp = 0.24. We 
have chosen rj = 10~^° (filled circles) and rj = 10~^ (squares) for purposes of illustration. 



expected behavior in the vicinity of the standard model, but we see that the curves turn 
over, and at fixed a there is a larger value of Hi for which the helium abundance is the 
same. This feature has been previously noted in the context of exploring the effects of 
very large numbers of light neutrino species (and thus Hi) ||2^ and in the explorations of 
BBN constraints on scalar-tensor theories of gravity |]2B|. For very large Hi, the freeze- 
out temperature is sufficiently high that n/p(Tf) asymptotes to unity; AteBN continues 
to decrease as Hi increases, leading to an ever-larger neutron fraction at nucleosynthesis. 
However, i^BBN also increases and eventually the expansion rate at nucleosynthesis is so 
large that there is not enough time to efficiently turn the neutrons into helium. Thus, 
given the success of standard BBN, there will be an additional larger value of Hi (at the 
standard value a = 2) for which the helium abundance is the same. This argument can be 
generalized not just to other values of a but also to the other light element abundances. 
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One expects (at fixed a) to get two values of i^BBN which produce the same abundance. 



4 Deuterium and Lithium constraints 

The deuterium abundance provides an important additional constraint on the allowed pa- 
rameter space. In the standard picture deuterium is a by-product of neutrons getting burnt 
to helium. Since the abundance of deuterium is set by its destruction rate and the time 
available for this destruction, it is extremely sensitive to the expansion rate at temperatures 
close to Tbbn- Around the standard model, D/H increases with Hi for a fixed a, since the 
time available for nucleosynthesis decreases. If one fixes Hi, the decrease in post-BBN 
(T < Tbbn) expansion rate causes D/H to decrease with increasing a. Once again we see 
that there is a trade-off to be made between a and Hi, and hence it is possible to increase 
a and Hi simultaneously such that we end up with the same deuterium abundance. This is 
borne out by the contours of constant deuterium plotted in Figure ^ Note that one expects 
exactly the opposite behavior with respect to changes in a (with Hi fixed), and Hi (with 
a fixed) for the upper branch of contour plot. Essentially, the abundances in the upper 
branch of the contour plot increase if the time available for nucleosynthesis goes up. 
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Figure 2: Contours of constant deuterium, D/H= 3 x 10 ^, at rj = 10 (filled circles) and 
Tj = 10~^ (squares). 
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From Figure H it is clear that the lower branch of the contour plot turns over at large a. 
This can be traced to the appearance of a new production channel. As a increases without 
a corresponding exponential increase in i^i, the time available for nucleosynthesis becomes 
very large, and consequently deuterium is almost completely destroyed. In the lower branch 
this implies that beyond a certain value of a it would not be possible to produce a reasonable 
D/H abundance. However, if the age of the Universe gets large enough it becomes possible 
to produce deuterium by pp fusion, in a manner similar to what happens in the sun. This 
pp burning explains the turn over towards smaller values of iJi in the lower branch of the 
deuterium contour. For this solution to work, the age during nucleosynthesis has to be of 
the order of billions of years. We, therefore, consider this solution highly unlikely and do 
not include it in our analysis. 

For our parameter space around the standard model we expect the same behavior with 
respect to ri as in standard BBN. All other things being the same, increasing r/ decreases 
deuterium abundance while increasing Yp due to more efficient nucleosynthesis. Thus at 
fixed a if ?7 is increased, then one needs to decrease ifi to lower the helium abundance back 
to where it was. In the case of deuterium, if 77 is increased at fixed a, then one needs to 
increase Hi in order to increase D/H back to where it was. 

Note that for the upper branch one would have to increase Hi at fixed a to compensate 
for an increase in 77 for hoth helium and deuterium, since both Yp and D/H increase with r\ 
in the upper branch. 

We can infer the primordial abundance of ^Li very accurately from poor-metallicity stars 



20| . Therefore, it provides an additional constraint on our parameter space. For most of 
the parameter space, getting helium right automatically ensures that lithium is near the 
observed abundance. Even so, with the very tight bounds on ^Li/H that we impose, lithium 
does provide strong constraints. Figure |^ shows contours of constant lithium in the Hi-a 
plane. Its behavior with respect to a change in rj follows that of helium. 

The ''Li iso-abundance contour has two interesting features not present in the helium 
contour, which can be traced to the fact that ''Li has two distinct production channels. 
^Li can be produced through the reaction ^He('^H,7) ''Li, and through ^He(^He,7)^Be with 
the subsequent beta-decay of ''Be. Fixing our attention on the lower branch of the lithium 
contour, as we move from low to high a, the dominant production channel changes from 
the direct ''Li channel to the indirect ''Be one. The kink in the lower branch of the contour 
is at values of a where this transition takes place. At large a and Hi, there is a small 
closed contour below the upper branch. For a fixed a in the region spanning the width of 
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Figure 3: Contours of constant lithium, ^Li/H= 3 x 10 ^'^j at 77 = 10 ^'^ (filled circles) and 
Tj = 10~^ (squares). 

the small closed contour, there are three values of Hi which produce the required lithium 
abundance. The solution with the lowest value of Hi produces ''Li through the indirect ''Be 
channel. In the upper branch of the contour ^Li production is through the direct channel, 
for all values of a. A larger Hi means shorter time for nucleosynthesis and since is easier 
to burn than ^He, the direct channel dominates. 



5 Combined constraints 

Putting together all of our constraints yields the allowed region shown in Figure ^. We note 
that the allowed region is essentially one dimensional, and is characterized by an almost 
linear relationship between log(ifi) and a. This tells us that there must exist a temperature 
Tc at which the expansion rate He is approximately fixed for all the models in the allowed 
range. In other words, we must have 

Hi = He {TJ MeV)" ; H^ = (0.039 ± 0.013) /sec at = 0.2 MeV . (9) 

The range of He quoted in Eq. ^ corresponds to our conservative limits on the baryon 
to photon ratio. For the case with the tighter bound on baryon density, we obtain H^ = 
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(0.03±0.004)/sec. The value of Tc reflects of the physics involved. ^He is primarily sensitive 
to Hi (as well as to a), while the other elements (which can be viewed as by-products of the 
neutron to helium burning process) are mainly sensitive to -f^BBN- Thus, not surprisingly, 
we have 1 MeV > Tc > Tbbn- 
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Figure 4: Allowed values of a and Hi which lead to acceptable light element abundances. 
Upper panel shows the allowed parameter space when we constrain rjiQ to be between 1 and 
10. The lower panel shows the allowed parameter space for 0.5 < rjiQ < 50. 

The coincidence that is so close to the annihilation temperature (~ me/3 = 
0.1 7 MeV) implies that the expansion rate of the Universe during annihilation must 
have been close to the standard BBN value. This in turn implies that the cosmic mi- 
crowave background will retain a black body spectrum through the annihilation epoch (as 
in standard BBN) in all the viable non-standard expansion histories we have identified. 

The effect of a change in the assumed range of the baryon density on the bounds on a 
is visible in Figure ^ Decreasing the lower bound on rj increases the upper bound on a, 
while allowing for larger baryon densities implies a smaller lower bound on a. The reason 
behind such behavior is not hard to understand. Given that the expansion rate at 0.2 MeV 
is approximately fixed, a viable expansion history with larger a has lower i^BBN- A lower 
i^BBN (i-e., more time for nucleosynthesis) at fixed baryon density implies a lower deuterium 
abundance. Thus one requires a smaller baryon density to get deuterium abundance right 
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in the limit of large a. Similarly, decreasing the lower bound on a requires larger baryon 
densities. 

As discussed before, we do not really have an observational upper bound on the baryon 
density of the Universe in the absence of a theory of gravity. The lower bound is on a more 
sound footing because we can count the visible baryons. Thus it seems pertinent to ask what 
happens if one allows for larger cosmological baryon densities. From our arguments above 
we would expect viable models with smaller values of a. However, increasing r/io above 50 
(panel 2 of Figure ^) has no effect on the allowed range of a. The reason behind this result 
is best understood through the following qualitative picture. At smaller values of a, and 
therefore larger i^BBN? one has a higher deuterium abundance if the baryon density is fixed. 
To lower the deuterium abundance to required levels one has to increase the baryon density. 
However, a larger baryon density (keeping a fixed) will lead to a larger helium fraction. 
Hence below a certain value of a (which depends on the constraints on the abundances) 
it is not possible to simultaneously obtain the correct amounts of deuterium and helium. 
This is also the reason why the allowed contour in the lower panel of Figure § does not 
appear to close. 



6 Neutrino asymmetry 

One of parameters in cosmology on which we have little handle is the neutrino asymmetry, 
i.e., the excess of neutrinos over anti- neutrinos or vice versa. This excess or deficit can 
be quantified by the neutrino chemical potential, /i, which enters into the Fermi-Dirac 
distribution function, given by 1/[1 + exp(p/T + yu/T)]. The chemical potential of anti- 
neutrinos is —fi. The ratio /x/T is an invariant, as long as there are no entropy-changing 
processes which involve neutrinos. 

There are no compelling experimental or theoretical reasons to expect the neutrino 
(lepton) asymmetry (oc (vr^ + fi"^ /T'^)fi/T) to be orders of magnitude different from the 
baryon asymmetry. In fact, one might very well argue that it would only be natural for 
the baryon and lepton asymmetries to be comparable if they were formed by the same 
processes. However, we have no concrete experimental hints to back up our theoretical 
prejudice. Moreover, in our framework, the value of fi/T is virtually unconstrained (but for 



BBN) since the standard cosmological bounds p3 do not apply. Note that in order to probe 



yu/T ~ 1 regime, one would require experiments which are capable of measuring neutrino 
energies down to 10"'^ eV. As we shall shortly see, BBN even without the assumption of 
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standard gravity provides us with constraints on the value of the electron neutrino chemical 
potential in the |yu/T| = 0{1) regime. 

Changing the neutrino chemical potential affects the neutrino number density. In stan- 
dard BBN, a change in the neutrino number density affects the expansion rate and hence 
the abundances. Clearly, this effect is non-existent in our study. The electron neutrino 
chemical potential also affects the n ^ p reaction rates. This is the effect we shall be con- 
cerned with in this section. From here on, /x/T will implicitly refer to the electron neutrino 
chemical potential. 

Detailed calculations of the effect of non-negligible fi/T (i.e., \fi/T\ > 1) on BBN have 
been performed p8[. The salient points are (1) the equilibrium value of n/p is modified by 
a factor of exp(— /x/T) from that when = 0, (2) the neutrino decoupling and n-p freeze- 
out temperatures are changed. The first point clearly implies the constraint, /x/T < 2, if a 
helium-4 mass fraction of 24% is to be synthesized. In Figure ^ we plot the allowed contours 
in a - Hi plane. The contours are labelled by the value of fi/T. We have restricted Yp to 
being 24.4%. Allowing for our previously adopted range in Yp will not change the result 
much as is apparent from the fact that the allowed contours in Figure | are practically 
1-dimensional. We note that if a is set equal to 2, then we are back to standard cosmology; 
Figure |^ then provides us with the constraint on the number of relativistic degrees of 
freedom (oc iJi) [^ . 

The principal feature in Figure ^ is the trend of the required value of Hi being larger for 
larger fi/T. At larger fi/T, the equilibrium value of n/p is smaller, and hence to compensate 
for that tBBN must be smaller or Hi larger. For fi/T < —10 no solutions were found, which 
owes to the fact that at the small values of Hi required to produce the right helium-4, ^bbn 
gets so large that one cannot synthesize acceptable amounts of D/H. 

7 Discussion 

We have shown that the viable (i.e., satisfying BBN abundance constraints) expansion 
histories can be well approximated by the range of allowed values of a, for a fixed value 
of fi/T. Lets concentrate on the more "natural" \fi/T\ << 1 (equivalent to /i = 0) case. 
Given the allowed range of a one can generate the allowed range of expansion rates at any 
temperature from about an MeV to 50 keV by using 

H{T) = H, {T/T,T . (10) 

We have plotted the history of these viable non-standard Universes in Figure ^. 
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Figure 5: Allowed values of a and for Yp = 0.244 and 0.5 < r/io < 50. The contours are 
labelled by the value of yu/T. No solutions consistent with our constraints were obtained 
for fi/T > 2 and fi/T < -10. 



The most interesting result of our study is that there is a range of expansion histories 
that are compatible with the observed light-element abundances. There is therefore room for 
substantial deviation from the standard cosmological model at early times, while remaining 
consistent with empirical evidence. On the other hand, it is encouraging to note that our 
allowed region (essentially one-dimensional) is only a small volume of the entire parameter 
space. In this sense it would be unlikely to find that any particular model was both very 
different from the standard picture, and consistent with the data. 

We have also considered the effect of a large neutrino chemical potential. Independent 
of general relativity, BBN constrains the value of the electron neutrino chemical potential 
to— 10<;u/T<2. Of course, this increases the allowed range of Hi considerably, as can 
be gauged from Figure ||. 



Throughout this discussion we have assumed tight but reasonable bounds |gO[ on the 
light element abundances. Our quantitative results are sensitive to the assumed ranges. As 
observations improve further, it will be possible to precisely pin down the expansion history 
of the Universe during the epoch of nucleosynthesis. 
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-10 -9.5 -9 -8.5 

log(a/ao) 

Figure 6: The expansion rate if as a function of scale factor a, for the family of allowed mod- 
els in our parameterization (corresponding to our conservative limits on the baryon/photon 
ratio and negligible electron neutrino asymmetry). The dashed curve represents the stan- 
dard model, and the solid curves are different allowed histories. It is clear that there is a 
value of the scale factor Cc/ao ^ 8.5 x 10~^° at which the allowed range in H is minimized; 
this corresponds to a temperature Tc ~ 0.2 MeV. 
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